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Abstract 

We present a model for the scaling laws of the electric field noise spectral density as a function of 
the distance, d, above a conducting surface. Our analytical approach models the patch potentials 
by introducing a correlation length, £, of the electric potential on the surface. The predicted scaling 
laws are in excellent agreement with two different classes of experiments (cold trapped ions and 
cantilevers), that span at least four orders of magnitude of d. According to this model, heating 
rate in miniature ion traps could be greatly reduced by proper material engineering. 

PACS numbers: 37.10.Ty,34.35.+a,72.70.+m,37.10.Rs 
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Device miniaturization is a challenge that raises new issues because the scaling laws valid 
in the macroscopic range might fail, for instance due to the emergence of a new characteristic 
length. Even the simple case of the electric field in the vicinity of a conductor surface can 
exhibit anomalous behavior caused by small inhomogeneities of the electric potential on 
the surface. These field fluctuations are crucial in the studies of short distance phenomena 
such as the measurement of the Casimir-Polder force [1], studies of non contact friction[2, 3], 
gravitational forces [4] and contact potentials [5, 6, 7]. 

In a different context, recent success in quantum information experiment with trapped 
ions (see Ref [8] for a review) motivated the fabrication of micro-traps in order to fulfill the 
scalability requirement of a quantum computer [9]. In such devices a set of micro-fabricated 
conducting electrodes generates an oscillating electric field that traps laser-cooled ions in 
a harmonic potential well, at a distance, d, of the surface. In this situation the presence 
of a fluctuating electric field affects the ion motion inducing a heating, that is usually 
characterized in term of quanta of vibration gained per unit time. This heating fixes a limit 
on the achievable fidelity of ion based quantum gates[10]. One might try to account for 
this heating rate by considering typical electric noise sources in conductors, among which 
the most likely is Johnson noise. However, measured heating rates are orders of magnitude 
larger than the expected contribution of the Johnson noise. Moreover, Johnson noise would 
induce a heating rate that scales as d~ 2 , whereas the observed one is consistent with a d~ A 
scalingfll], as would be expected from a random distribution of charges, leading to the 
notion of patch potentials[12]. 

Recent experiments[13] suggest that indeed the surface quality plays a dominant role in 
this anomalous heating. The observed scaling of the field noise with temperature points out 
thermally activated phenomena (surface defaults, charge traps,...). These effects are easy to 
probe in the static limit where direct observation of the surface is possible using an atomic 
force microscope, but much more difficult to observe directly at higher frequencies, which 
are relevant to ion trapping. However the measured noise can be compared by assuming a 
variation as c^ -1 , experimentally verified at room temperature [13]. 

The question of finding the electric field fluctuations near an infinite conductor filling half 
the space has already been extensively treated by considering thermal fluctuations carried by 
uncorrelated punctual sources (Johnson noise) [12, 14, 15]. In these models the characteristic 
length is given by the skin depth of the material, S, at the considered frequency. The scaling 
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of the electric field noise density is then expected to change from (d/5)~ 3 for d <C 5 to 
(d/5)~ 2 for d 5 where the typical value of S for gold electrodes in the MHz range is tens 
of microns. However this behavior has not been observed in ion traps experiments where 
a d~ A scaling is reported for d varying from 75 \im to 1 mm[16]. In a recent work [17], this 
type of model as been improved to take into account effects of charge diffusion, where the 
characteristic length, 5, is related to the mean free path of charges on the surface. The 
predicted scaling is consistent with (d/S)~ 3 for d <C 5 to (d/5)~ 4 for d ^> 5. 

In this work we present a model where we take explicitly into account the spatial 
dependence of the electric field noise density above a surface, introducing a characteristic 
length not yet considered in previous models: the correlation length, (, of the noisy po- 
tential on the surface. In that picture, the noise arise from finite size sources distributed 
randomly on the surface. We show that this simple model accounts for the behavior of elec- 
tric field noise scaling on the whole range covered by both ion heating measurements and 
cantilever-based measurements. We propose that this correlation length could be related 
to the characteristic size of the patches, opening a way to the control of noise intensity by 
engineering the material properties. In a recent work, an analytical solution of the Laplace 
equation for the special case of a planar ion trap was found [18]. In the present letter we 
extend this approach and obtain an analytical expression for the scaling law of the electric 
field noise density Se(&, d) at a distance d of the surface. The first step is to solve the 
Laplace equation for the potential </>(x,y,z,t), A<f>(x,y, z,t) = in half the space (y > 0) 
with boundary condition <j)(x, 0, z, t) = 4>o{x, z, t) and vanishing potential for y — > oo. Under 
these assumptions, the potential reads: 



where we took the Fourier transform of the Laplace equation, introducing k x and k z . Then 
carrying out the integration on k x and k z , the resulting potential above the plane is: 



where we introduced K(x, y, z) = ^ x2+y 2 +z 2ys/2 , the kernel of the Laplace problem with these 
specific boundary conditions. As in ref [18], equation (2) allows us to compute efficiently 
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the potential associated to a given electrode geometry. More generally equation (2) allows 
us to compute the potential created by any disordered boundary condition <f) (x,z,t). 

In the following we will assume that small patch potentials distributed over the plane 
(x, z) create a disordered electric potential on the surface. Let {Cj} be the area of these 
patches and Vi(t) the time-dependent electric noise on the patch Cj. With these nota- 
tions one have: (po(x, z, t) = 'Yl li Vi{t)xc i {x, z) where xCi is the characteristic function of Cj 
(.XCi( x i z ) — 1> z) G Ci and Xd( x , z ) — 0, V(x, z) Cj). From equation (2) we can 
compute the electric field temporal correlation function at a distance d: 

s E (r,d) = ^^IMM 

i 3 

x / du'dv' XCl («', v') [VK] (x - u, d,z- v') 



x J dudvxci{ u -, v ) [VK] ( x — M ? di z ~ v ) (3) 

where [VK] is the gradient of the scalar field K(x,y,z). The horizontal line in equation 
(3) means that we average over many configurations of Vj,. We suppose that the noise on 
two distinct patches originates from independent random processes with the same temporal 
correlation function R{t): 



K(fM(f + r) = ayi2(T) (4) 
Equation (3) can then be rewritten: 

S E ( T , d) = / dudvdu'dv' ^ Xc t (u, v)xc\ («', «') 

x [VK] (x-u : d,z- v). [VK] (x - u', d,z- v') (5) 

Since the set {Ci}i is disordered one can write: ^ XCiiu, v )xc t (V, v') = 
N (XCi{ u ,v)xCi{ u ' ' i v ')) 1 where the brakets denote an average on the configuration of {C{\ 
and N is the total number of patches. Introducing the spatial correlation function of the 
patches, C^(u — u',v — v') = (xCi(u, v )xcA u 'i v '))i where ( is the (finite) correlation length, 
we take the temporal Fourier transform of Se{t, d): 



2 

i(k x u+k z v) 



(6) 
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where Sv(cu), defined as the Fourier transform of -R(r), is the potential noise spectral density 
on the surface and S^(k x ,k z ) is the two dimensional Fourier transform of C^(x, z). After 
some calculations detailed in the Appendix, one obtains: 

S E (co, d) = NS n V ;> [ dkdOSAk cos6, k sm6) k 3 e~ 2dk (7) 

21T Z J 

To simplify further this equation we need an explicit form of S$. In what follows we assume 
an exponential behavior for the spatial autocorrelation function of the patches: 

C c (x,z) = e - v/ * 5 + I X (8) 

This correlation function arises from a Poisson Voronoi tessellation model of the polycrys- 
talline structure of metals[19]. More generally, following the approach of Ref [20], it can be 
demonstrated that such an exponential behavior arises from a collection of random potential 
patches on a surface. 

We can then compute the Fourier transform explicitly: S^(kcos9,ksin9) = n+pfsw^ ; 
finally giving: 

N( 2 f°° k 3 e~ 2k 

S E (u : d) = 2—S v (u)J o dk {d2 + ek2)3/2 (9) 

Equation (9) clearly identifies the separate contribution of the spatial and temporal com- 
ponents of the potential noise on the plane to the electric field noise density at a distance 
d of the surface. Let us point out two important limits: in the case d 3> (, one finds 
Se{oj,<J) ~ — and in the case d <C C> one finds ^(w, d) ~ ■ Between these two 
simple asymptotic behaviors a smooth transition occurs as shown in Fig. 1. Let us note that 
the choice of any correlation function in (8) would have led to the same long range limit 
(d 3> C), whereas this choice is crucial for the short range behavior. 

Let us now analyze how equation (9) applies to both the case of ion traps and cantilevers 
based experiments. In planar ion traps, charged particles are trapped above the surface in 
a tight pseudo-potential and are usually cooled down to the fundamental vibrational state. 
The presence of a fluctuating electric field leads to the heating of the ion with a heating rate 
r = j^-t- Se(uj, d), where u is the vibrational frequency of the ion (typically a few MHz), 
m is the mass, e the charge and d the equilibrium position of the ion[12]. As described 
in Refs[13, 21], the experimental measurement of ion heating rate gives thus access to the 
electric field noise density. Due to electrode configuration in planar ion traps, the ions are 
trapped at a distance d proportional to the size of the electrodes: the field can thus be 
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FIG. 1: (Color online) continuous line: normalized field noise density Se(w, d) (in units of NS v^) ^ 
as d varies for a fixed £. For d <§C ( one has SE(to,d) ~ d^ 1 (dashed line). For rf > ( one has 
SE{oJ,d) ~ d~ 4 (dotted line). 

probed in different ranges (typically [75 — 150] /xm) using different traps[13]. Although the 
surface quality (and thus () might depend on the fabrication process, a d~ A dependence of 
the heating rate (or field noise density) has been observed. These results are in agreement 
with the limit d 3> C of our model. 

In cantilever based electric field noise measurement, the potential above a surface can 
be probed on a typical distance range from 10 nm to 400 nm, yet on a very different scale 
from ion traps. The cantilever is a device with resonant mechanical oscillation frequencies. 

2 

Its movement is damped by coupling to stray electric fields, with a rate T = j^;Se(uj, d), 
where |^ is the frequency of the cantilever (typically a few kHz) and q = CV is the induced 
charge, equal to the tip-sample capacitance, C, times the potential bias, V[3]. Measuring the 
cantilever oscillations damping rate using optical interferometry gives access to the electric 
field noise density[2, 3]. These measurements, in the range d ~ [10 — 100] nm, give a d~ l 
scaling of the field noise density, consistent with the limit d C, of our model. 
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As mentioned above, the frequencies probed in ions traps and cantilever based mea- 
surement differ by several orders of magnitude, but the measured noise can be scaled to the 
same frequency, assuming a variation as u" 1 . 

Table I summarizes electric field noise density measurements above gold surfaces, as 
reported in several experimental works. 



Ref. 


d 


w/2tt 




S E {uo,d) 




lim 


MHz 


V 2 m~ 2 Hz~ 1 


[2] 


0.02 


4 x 10" 3 


4 


1.6 x 10~ 2 


[22] 


40 


3 


9 x 10~ 12 


2.7 x 10" 11 


[13] 


75 


1 


[0.3 - 3] x 10- 11 


[0.3 - 3] x 10- 11 


[12] 


140 


10 


5 x 10" 12 


5 x 1Q- 11 



TABLE I: Experimental room temperature SE(uJ,d) values taken from cited references. Above the 
horizontal line can be found a value measured using a cantilever. Below this line lie values measured 
in ion traps experiment. The last column give the rescaled value of <Se(u;o, c?) for ojq/2tt = 1 MHz. 
Data from reference [13] have been obtained with the same trap after successive cleaning procedures 
and thermal cycling. 

Let us compare our model to the rescaled measured values, S#(<i), of Table I that we 
report in Fig. 2. In order to plot on the same graph the field noise density given by equation 
(9), we need numerical values for NSy(uo) and £■ As the cantilever data are consistent 
with the short range limit of our model, Se{&o, d) = NSv{uJo)/{(d) we can obtain a value 
for NSv(ujq) = 3.2 x 10~ 10 Co, where we introduced £o the characteristic length of this 
sample. Under the simple assumption that NSv(uj ) is not sample dependent, we can plot 
the curves corresponding to different values of £/Co (dashed lines in Fig. 2). We find that 
0-6£o < C < 4.5(o covers the range of ion traps data. As reported in Refs [5, 7], the typical 
size of static patch potentials on gold surfaces is (q = 1 fim in excellent agreement with 
the assumption (q 3> 20 nm. We interpret then the residual spread of the experimental 
data in terms of different correlation lengths (in the range [0.6,4.5] fim) associated to the 
sample surface quality, itself highly dependent on the fabrication technique. The role of 
surface quality has been noted in ref [13], where successive cleaning procedures decreased 
the measured electric field noise density (points Q m Fig- 2). Thus our work shows that in 
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FIG. 2: Plot of experimental data from Refs. [2] (■), [22] (0), [13] (O) and [12] (□). The lines 
correspond to the value of Se(loo, d) predicted by our model for three values of £, Q = 0.65Co (solid 
line), C = l-6(b (dashed line), £ = 4.6Co (dotted line). 

order to compare measured field noise on very different scales, one has to take into account 
the correlation length of the material and not assume a d~ A scaling of the noise, otherwise 
leading to an over-estimation of the noise at smaller distances. For example, as shown in 
Fig. 1, our model predicts that for a sample with £ ~ 1 fim, the noise at d = 100 nm is 10 9 
times greater than noise at d = 100 fim (at the same frequency), three orders of magnitude 
lower than expected with a d~ A scaling. 

In conclusion we developed a model for the electric field noise density based on an ana- 
lytical approach, assuming a finite size of the potentials patches. Electrical noise field density 
scaling laws predicted by this model are in very good agreement with experimental results, 
both in the long range regime (hundred of microns), probed with ion traps experiments and 
in the short range regime (tens of nanometers) probed with micro-fabricated cantilevers. 
Let us point out that the pessimistic d~ 4 scaling law for heating rates observed in trapped 
ion experiments is over estimated in the short range regime. Moreover even though the 
noise considered in such devices lies within the MHz range, we find a characteristic length 
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compatible with measurement of static potential patches sizes. This opens a new possibility 
to improve the performances of surface micro-traps based on the analysis of electric static 
noise obtained with atomic force microscopes. 

We thank E. Boulat for fruitful discussions. This work was supported by the French Na- 
tional Research Agency (ANR) Project No. ANR-JC0561454. R.D. gratefully acknowledges 
the funding from the Delegation Generale de l'Armement (DGA). 
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Noticing that K(d x x,d x y,d x z) = d 2 K(x, y, z) we change variables in equation (6) 
so that: 

Ssfad) = [ kdkd8S c (k cos 6, k sin 9) 



x 



16vr 4 ci 2 



ududcf) \VK] (u cos <f),l,u sin <p) e 



-idku cos(0— 8) I 



where u is now dimensionless. Since: 



\VK] (u cos 4>,l,u sin < 



Se{u, d) may now be rewritten: 



/ 



(l +w 2)5/2 



— 3u cos 
u 2 -2 
—3u sin 



A 



S E{u,d) = N A S j}^ I kdkd9S c (k cos 9, k sin i 



4vr 2 rf 2 



du , \^ /2 

(1 + M 2 ) 5 / 2 



/ \ 2 

' -3m cos 9 Ji (dku) » 

(u 2 -2)J (dku) 
^ —3iu sin 9 J\ (dku) j 



(A.IO) 



J n (x) being the n-th order Bessel function of the first kind. The integral on u can be 
analytically calculated leading to equation (7). 
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